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Abstract. Electrical impedance tomography (EIT) is an imaging modality 
in which the conductivity distribution inside a target is reconstructed based 
on voltage measurements from the surface of the target. Reconstructing the 
conductivity distribution is known to be an ill-posed inverse problem, the solutions 
of which are highly intolerant to modelling errors. In order to achieve sufhcient 
accuracy, very dense meshes are usually needed in a finite element approximation 
of the EIT forward model. This leads to very high-dimensional problems and often 
unacceptably tedious computations for real-time applications. In this paper, the 
model reduction in EIT is considered within the Bayesian inversion framework. 
We construct the reduced-order model by proper orthogonal decompositions 
(POD) of the electrical conductivity and the potential distributions. The 
associated POD modes are computed based on a priori information on the 
conductivity. The feasibility of the reduced-order model is tested both numerically 
and with experimental data. In the selected test cases, the proposed model 
reduction approach speeds up the computation by more than two orders of 
magnitude in comparison with the conventional EIT reconstruction, without 
decreasing the quality of the reconstructed images significantly. 
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1. INTRODUCTION 

In electrical impedance tomography (EIT), a set of alternating currents are injected to 
the target through an array of boundary electrodes, and the resulting voltages on the 
electrodes are measured. Based on these measurements, the internal conductivity 
distribution of the target is reconstructed. The applications of EIT include e.g. 
industrial process monitoring, [121 [Ml 133 US] , control [HI and design [JS], medical 
imaging [TOl [H] and geophysical exploration [H [TTl [35]. Advantages of EIT are the 
very good temporal resolution (up to 1000 frames per second |48j). the non- invasive 
and radiation-free nature, and affordable measuring devices. A drawback of EIT is 
often the relatively low spatial resolution, which is partly due to the diffusive nature 
of the modality and partly a result of simplifications that are typically made in image 
reconstruction in order to compute the reconstructions in the required time frame. 

The image reconstruction problem in EIT is known to be an ill-posed inverse 
problem, i.e. the solutions are very sensitive to measurement noise and modelling 
errors. Therefore, an accurate model for the measurements - i.e. the forward model - 
is needed in order to obtain feasible reconstructions. The most accurate forward model 
for a physically realizable EIT measurements is referred to as the complete electrode 
model (CEM) [7]. The CEM consists of an eUiptic partial differential equation (PDE) 
and associated boundary conditions. The CEM is usually approximated with the 
finite element method (FEM) where finite dimensional approximations to the electrical 
conductivity and potential are written using locally supported piecewise polynomial 
basis functions. Typically, dense finite element (EE) meshes are needed in order 
to model the measurements with sufficient accuracy. Therefore, the use of locally 
supported EE-bases often leads to a high dimensional problem which may require a 
large amount of memory and time to solve. The image reconstruction problem of EIT 
is a non-linear inverse problem, and in optimization based reconstruction methods, 
the conductivity estimate is computed iteratively. In the case of a conventional EE 
approximation of CEM, the large dimensional forward problem needs to be solved 
in each iteration step. However, in many industrial and medical applications online 
reconstructions are needed, and the allowable time for the computations is short - 
sometimes even in order of milliseconds. Usually, the problem of high computational 
demand is overcome by using coarse FE-meshes j2H 133) and/or global linearization of 
the observation model [45 [137] or even highly simplified back projection reconstruction 
methods [3l[40]. Such choices, however, often lead to severely biased reconstructions. 
It is thus of great importance to develop new methods for model reduction in EIT, 
so that the computation time could be decreased without ruining the quality of the 
reconstructions. 

The proper orthogonal decomposition (POD) has been widely used in different 
types of data and model reduction problems. In data reduction, data vectors are 
represented in a low dimensional space spanned by POD modes which are eigenvectors 
of the covariance matrix of the data set. The data reduction applications of the 
POD include e.g. image processing [51], face recognition [9l [24l [41], and data 
clustering [5^. Further, the POD has been applied to the reduction of PDE models 
[6l m [HI [311 [39] . Ii^ these papers, the basis functions for the EE approximation of the 
PDE were constructed by using the POD. The POD modes were computed based on 
the conventional EE approximation of the PDE with some fixed parameters. 

The POD and other projection-based model reduction methods have also been 
applied to inverse problems spanned by PDEs. Firstly, the POD has been applied 
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to constraining the solution of the inverse problem (i.e. the unknown parameter 
distribution) to a desired subspace which is selected based on the known properties 
of the target. In EIT, the POD modes for the internal conductivity distribution have 
been selected based on anatomical information in imaging of thorax |44| . Respectively, 
in [TJ [5] the POD was used to construct a reduced-order basis for the unknown 
material parameters in other electromagnetic inverse problems. Secondly, the POD 
has been used for reducing the number of basis functions for the state of the system 
(i.e. the solution of the PDE). For example, in [17] the temperature distribution 
was represented in a reduced-order POD basis in an inverse heat conduction problem. 
Here, the inverse problem was to estimate a scalar Robin coefficient using temperature 
measurements from the boundary of the domain. An ensemble of temperatures for 
constructing the POD basis was obtained by solving the heat equation corresponding 
to an ad hoc selected set of Robin coefficients. For other studies where the POD modes 
were used for constructing the FE basis for the solution of a PDE, see [SJ |H1 HH 1311 ISS] ■ 
Note that in all these papers, the number of unknown parameters in the inverse 
problem was small, and the POD was applied only in the state of the system. Recently, 
Lieberman et. al. |28| proposed an approach to the model reduction for high- 
dimensional statistical inverse problems. They applied a heuristic optimization-based 
greedy sampling for simultanously constructing the bases for both the unknown high 
dimensional parameter distribution and the state of the system. The feasibility of the 
proposed approach was demonstrated with an example of a statistical inverse problem 
of groundwater flow. By using the reduced model, they computed an estimate for the 
posterior distribution of the hydraulic conductivity parameters by a Markov Chain 
Monte Carlo (MCMC) method. 

In this paper, we propose a POD-based reduced-order model for EIT. As in 
[28j . we represent both the parameters (the electrical conductivity distribution) and 
the state of the system (the electrical potential distribution) in reduced bases. We 
consider the model reduction within the Bayesian framework, and write an explicit 
statistical prior model for the unknown conductivity distribution. Unlike in [28], we 
construct the reduced order basis for the conductivity based on the prior distribution. 
Moreover, we compute the realizations of the potential field using samples drawn from 
the prior distribution of the conductivity, and use the ensemble of potential fields to 
construct the POD basis for the potential. Finally, we complete the reduced model by 
approximating the errors related to the reduced-order approximations by an auxiliary 
additive noise process. In the last step, the so-called approximation error approach 
[H m] is adopted. 

The rest of the paper is organized as follows. In Section [21 the EIT imaging 
is briefly reviewed; the discussion is limited to the computational aspects of the 
forward model and the Bayesian inversion. The POD based model reduction for 
EIT is proposed in Section [3l Further, in Section HI the reduced-order reconstruction 
is evaluated both with simulated and experimental data. Finally, the conclusions are 
drawn in Section [5l 

2. ELECTRICAL IMPEDANCE TOMOGRAPHY 

In this section, computational aspects of EIT imaging are considered in Bayesian 
(statistical) framework. The complete electrode model (CEM) and its FE 
approximation are briefly reviewed in Sections 12. II and r2. 2) respectively. The Bayesian 
inverse problem of EIT is discussed in Section 12.31 
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2.1. Complete electrode model 

In EIT, the most accurate measurement model is the CEM [7], which consists of the 
PDE 

v-(CTVit) = o, fen (1) 

and the fohowing boundary conditions: 
Qu 

u + zea—^Ui, fe e£, £ = 1,2, . . . ,i (2) 
on 

=0, fedn\utiei (3) 

on 

3n 

a—dS = Ii, e = l,2,...,L (4) 



d 



n 



where a ~ is the electrical conductivity, u ~ u{f) is the electric potential 

inside the target domain fJ, also referred to as the inner potential, f is the spatial 
coordinate and n is the unit outward normal vector. Contact impedances, electrode 
potentials, and injected currents corresponding to the electrodes eg, £ — 1,...,L 
are denoted by ze, Ui, and If, respectively. For further use, we denote the vectors 
z = [zi, . . . U = [Ui,...,Ul]^ and / = [/i, . . . , /l]'^. In addition to (HH!]), we 



write 



L L 

^/, = 0, ^[/, = 0. (5) 



i=i e=i 

where the former condition ensures that the charge conservation law is fulfilled and 
the latter one fixes the reference potential level. 

The forward problem of EIT is to solve the inner potential u = u{f) and the 
electrode potentials Ui, i = given the conductivity distribution a{f), the 

contact impedances zg and the injected currents Ig. The existence and uniqueness 
of the solution was proven in |42| . The measurements in EIT imaging consist of 
potential differences, i.e. voltages, between electrodes. The inverse problem in EIT 
is to reconstruct the conductivity distribution a{f) given the voltage measurements 
corresponding to various sets of electrode currents. 



2. 2. Finite element approximation of the CEM 
The variational form of ((THll) can be written as [42] 

L 

B{{u,U),{v,V)) = Y,IiVi , 'i{v,V)eH (6) 
(=1 

where H = H^{il) x R^, H^{Q,) is a Sobolev space and B : H x H ^ R is a bilinear 
form such that 

B{{u,U),iv,V)) = [ CTVit ■ Vvdx + V - / iu-Ui){v-Vi)dS{7) 

Jn Jei 

In this section, the FE approximation is reviewed briefiy. For details, see [45]. Later, 
in Section[31 we point out how the reduced order approximations affect the FE scheme. 
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In a finite element implementation of the CEM, the finite dimensional 
approximations for the conductivity and the inner potential u are written as 

N M 

a^a'' = J2 a^cj,, (f) , w « = ^ /3,^, (f) . (8) 
1=1 ]=i 

Here = a^{r) denote the approximated conductivity and inner potential 
respectively, (j>i = (j)i{r), i = 1,2, ...,A'^ and ipj = ipj{r), j = 1,2,...,M are 

the basis functions, and a = [ai, • • • , ckA/]"^, /3 = [Pi, 132, ■ ■ ■ , f^Ai]^ are the 
corresponding coefficients. In the discussion below, we identify a and its finite 
dimensional representation a. Usually, the basis functions 0^ and -0^ are selected 
as piecewise linear or piecewise quadratic functions. For higher order polynomial 
approximations, see |3H]- The potentials U on the electrodes are written as 

L-l 
k=l 

where Uk are the basis functions chosen as ni = [1, — 1, 0, . . . , 0]""", n2 ~ 
[1,0,-1,...,0]'^, riL-i = [1,0,0, ...,-1]"^ to fulfin the latter condition in ©. 

Inserting the above approximations to the variational form ([6]) leads to the 
following matrix equation |45) 

A9 = f (10) 
where 9 = [/S^ , j^]'^ G r(m+l-i)xi^ ^ ^ [71, . . . , 7i_i]T and the vector / G 



L i)xi is defined as 
/ 



Om 



(11) 



Here, Om E R^-'^'^ is a vector of zeros and C = [ni, . . . , tt-l-i] G R^^'-^ ^\ 
Furthermore, matrix A = A{a, z) G R(A/+i-i)x(A/+L-i) ^^^^le form 



A = 



B + D EC 
C^E"^ C^FC 



(12) 



where 



^(*' j) = / 'l'^^^^ ■ VV-jdr!, 1 < 1,3 < M (13) 

^=l 

^(^'J) = E - / ^^^J^^' 1 < < M (14) 



E{i,j)^ -- [ ^,dS, l<i<M,l<j<L (15) 

Aw f 0, i^j 

F{hj) = E - / {^Mn,)tdS =\ , - , 1 < < L (16) 

where \ej\ is the measure of the electrode Cj, i.e. the length of the electrode in two 
dimensions (2D) and the area of the electrode in three dimensions (3D). The integrals 
are usually computed numerically for example with Gaussian quadratures. 

Formally, solution of the system pO|) is 6* = A~^f. In practice, however, 
the system of equations is solved, for example, with the LU-decomposition. As 
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noted above, the solution 9 consists of the parameter vectors /3 and 7 corresponding 
to approximations of internal potential field u and the electrode potentials U, 
respectively. This is the solution of the forward problem. It should be noted that 
if the basis functions ipj for the internal potential are locally supported, matrices 
B and D are sparse. Further, the demand for the high accuracy approximations 
often necessitates the use dense FE meshes which implies that B and D are high- 
dimensional. On the contrary, if the potential u is approximated with a low number 
of globally supported basis functions - choice made in Section [3]- the matrices B and 
D become low-dimensional and dense. 

Based on the finite element approximation described above, a set of voltages V 
between selected electrodes can be written in the form V = MCj, where 7 = 7(0-, z, I) 
is obtained from the solution of ([TU]), and is a difference matrix referred to as the 
measurement pattern. Further, since the dependence of the current pattern / and 7 
is linear, we can write 

V = R{<T,z)I (17) 

where R = R{<7, z) is referred to as the resistance matrix, and i?(cr, z)I = A4C"f. 
In EIT, the voltages V^^*^ are measured corresponding to several current patterns 
/^*\ i = 1, . . . , iViiij, where A'inj denotes the number of current patterns. If the 
measurement pattern A4 is the same for all current injections, equation (|17p yields 
= i?(cr,z)[/(i),...,/(^'"i']. For the notational convenience, we 
assume here that the contact impedances z are known; however, the extension to 
cases of unknown z is straightforward, see |37]. Then, by stacking all the voltage 
measurements in one column vector = [(y^^^)""", . . . , (l/(^inj)^TjT^ write the 

observation model of EIT 

V, = V(a) + e (18) 

where V(ct) = [(i?(tT, z)/(i))'^, . . . , (i?(cr, z)/(^'"i))'^]'^ and e is a stochastic term 
consisting of the measurement noise. Note that the observation model ([18]) is non- 
linear with respect to the conductivity distribution a. 



2.3. Bayesian inversion in EIT 

The solution of the statistical inverse problem in EIT is the posterior distribution, 
i.e. the conditional distribution of cr given the voltage measurements. We denote 
the associated posterior density by TT{a\V^^°^) where VT^^^** is the realization of the 
voltages Ve . Using the Bayes formula, the posterior density can be written in the form 

(-T/inoasI ■\ ( \ 

H^\Ve ) = (-^mcas) ^^i^e (19) 

where tt{V^^''^^\(t) is the likelihood density defined by the observation model ([T8| . and 
T^Vc{V^^'^^) — /""(o"! V^^°^^)d(j acts as a normalization factor. Further, na-icr) denotes 
the marginal density of ct, also referred to as the prior density because it includes the 
information we have on the conductivity before the measurements. 

A practical estimate for the conductivity is obtained by computing some point 
or spread estimate from the posterior distribution. One of the most commonly used 
point estimate is the maximum a posteriori (MAP) estimate 

(JMAP = arg max 7r(a| 1/'"'='^''). (20) 
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The problem of finding the MAP estimate requires solving an optimization problem. 
For example, if a and e are modeled as mutually uncorrelated Gaussian random 
variables cr ^ A/'(<t, Fg-) and e ~ A/'(e, F^), the MAP estimate gets the form 

auAP = argmin - a)f + - V(a) - e)\\^} (21) 

where Lo- and are the Cholesky factors of the inverted covariance matrices F~^ 
and F^^, respectively |20| . Due to non-linearity of the mapping V(cr), the solution 
of the minimization problem (j2ip is sought iteratively by some optimization method, 
such as the Gauss-Newton method. In each iteration step, the function V(cr) and the 
Jacobian matrix J = '^'^J'^^ are evaluated. For computation of the Jacobian matrix 
with the so-called adjoint method, see |19) . 

It is worth to notice that since the conductivity distribution cr is a random 
variable, the internal potential distribution u is also stochastic. This realization is the 
core of the next section in which both the conductivity and the potential distribution 
are represented in reduced-order bases selected with the POD. 



3. REDUCED-ORDER MODELLING IN EIT BASED ON PROPER 
ORTHOGONAL DECOMPOSITION 

In the POD, samples of (possibly) correlated random variables are converted into a 
set of uncorrelated random variables. In this procedure the original data is projected 
into a subspace spanned by orthogonal basis functions referred to as the POD modes. 
Depending on the application, POD is also known as principal component analysis 
(PGA) p^, Karhunen-Loeve decomposition [531 [551 [3D] , or Hotelling transform [TF| . 
The governing idea of POD is that in a case of highly correlated data, a low number 
of POD modes account for most of the data. For this reason, POD is often useful for 
data and model reduction problems. 

In the POD based reduced-order model the FE approximations cr'* and m'' (see 
equation ([5])) for the conductivity and potential are further approximated as 

a«a''«ao+f;afOD0POD(,-^ (22) 

i=l 
M 

^/3;°°^JO°(r) (23) 



where erg = E{cr} and uq = E{m} are the expectations of cr and u, respectively. 
Further, (pf^^ and ipf'-'^ are the basis functions of cr and u, and af^^ and /3j'^^ are 
the corresponding coefficients. The basis functions (^f-*^ and f/'f^^ are selected such 
that E{||a'* - ao - e£i «r°° '^f^^lP} and E{\\u'^ ~ u, ^ Y.f=i /^^^ V-r^H'}' ^re 
minimized over all N and M dimensional bases, respectively. Such a representation 
is referred to as the POD |18) . 

In practice, the POD bases {(/>f*^^} and {V'J''^^} of dimensions N and M 
are obtained as eigenvectors of conductivity and potential covariances T„ and F„ 
corresponding to N and M largest eigenvalues, respectively. Here, the covariance 
matrix T„ is determined by the prior model Tr^ (cr), and the covariance of the electric 
potential, F„, is approximated by the sample covariance 

1=1 
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where u^'^ is a sample from the potential distribution and = y ^J=i ^''^ 
the sample mean. The samples m^'^ are computed by solving the forward problem 
corresponding to samples {ct^*^}^_j^ drawn from the prior distribution of a. That 
is, corresponding to each conductivity sample cr^*), the CEM (HH?]) is approximated 
with the FEM as described in Section [521 Here, the conventional locally supported 
piecewise polynomial FE-bases {4>ii^} and {ipji^} for f and u are used. 

Denote the eigenvalues and eigenvectors of by Xf and vf , respectively, and 
those of r„ by and v^. The POD bases {«!>r°°(r)} and {ipj'^^ir}} are 



/POD 



(r) 



N 

E 

k=l 
M 

E 

k=l 



(25) 



(26) 



Due to the construction, we have for the reduced order conductivity parameters 

POD _ (^POD^ 



a 



(27) 
(28) 



r^POD =diag(A?,...,A^) 

where F^jPod denotes the covariance matrix of a^*-"^. 

The reduced-order observation model of EIT is constructed by inserting the 
approximations (|22p and ([23)) into the variational form ([6|). This yields, equivalently 
with ([To]) , a system 

^POD^POD ^ f (29) 
^,/3POD^[^POD^,_,^^PODlT 



Here 0^°° = 



^podT^^t 



A 



POD 



^POD 



^POD 
T 



M 

^POD^ 

C^FC 



where 



B 



POD 



JV 

E' 

fc=0 



POD rPOD 



L 

E 



< i,j < M 



^pod^pod^^^ 0<z,j<M 



POD/ 



1 



ihj) = - - / V'r°°d5, < z < A/, < J < L 



(30) 

(31) 
(32) 
(33) 
(34) 



and f,^,Zj,C and F are defined in Section 12.21 Here, for notational convenience, 
the expectations of a and u are denoted by ^q*-*^ and ip^'^^. Hence, 0Q'-*^(r) = 
(70, a^°° = 1, ?Ao°°(^) = "0 and 13^^^ = 1. When solving the system ([29]), the 
terms depending on (3q'^^ are moved to the right hand side of the equation. 

Because the POD bases used in the reduced-order approximation of u are globally 
supported, matrix A^'~'^ is dense. However, if the number of selected POD modes (A/) 
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is small, the dimension of A^'~'^ is low, making the solution of the reduced-order CEM 
computationally inexpensive. The same applies to computing the Jacobian matrix of 
the reduced-order observation mapping. 

We complete the reduced-order observation model by accounting for the error 
caused by the model reduction. Here, we adopt the approximation error method [H] 
which is based on statistical modelling of errors. Previously, the approximation error 
method has been used for recovering from errors caused e.g. by discretization |2HI33|. 
uncertainty of the geometry |13| . unknown boundary data |27| and unknown contact 
impedances |321 [55] . Here, we rewrite the observation equation ^TE\i in the form 

T/, =vPOD(aPOD)+e" (35) 

where V^'-'^ (a^'-'^ ) denotes the reduced-order counterpart of the mapping V{a) in 
psp. Further, the error term e" is of the form e" = e + e' where e is the measurement 
noise, and 

= V(a)-VPO°(aPO°) (36) 

is the error caused by the reduced-order approximations. To construct a model for the 
statistics of the error term e', we compute e'^^^ = V((7(*^) —V^'~'^{aJ^^^) corresponding 

T 

to all conductivity samples cr(*) in the sample set {cr^*-'}^ .j^. Here ceji^^, a sample of 
the reduced order representation of conductivity, is computed by projecting (T(,;) to the 

reduced order subspace spanned by {<P^^^}^^i- We approximate the expectation and 
the covariance matrix of e' by the sample mean and covariance, and write a Gaussian 
approximation for e', such that e' ^ 7V(e',re/). Further, if e' and the measurement 
noise e are mutually uncorrelated, we can write e" ^ A/'(e", T^n). where the expectation 
and the covariance of e" are of the form e" = e + e' and F^// = F^ -I- T^i , respectively. 

Using the reduced-order models, we write the MAP estimate corresponding to 
(|2ip in the form 

«MA? = arg min {||i,PODaP°°||2+ 

||^^„(^mcas _ ^POD(^POD) _ -/,)||2| (37) 

where we have noted that according to equation (|?7)) the expectation of the parameter 
vector Eja^*-"^} = [0, . . . ,0]'^. Further, Lq,pod and L^n are the Cholesky factors of 
the inverted covariance matrices F^poo and F^,^, respectively. The form (|57|) implies 
that e" and a^*^^ are mutually uncorrelated; this approximation has turned out to be 
adequate in many cases [21] . For modeling the cross-correlation of e" and 

aPOD^ see 

m- 

Finally, it should be noted that constructing matrix A^*-"^ in the reduced order 
EIT forward model is a time consuming task, especially because it necessitates 
the sample set {u^*)}. This set is constructed by solving the (original, high 
dimensional) forward problem of EIT T times, and the number of samples (T) 
must be set large in order to retrieve the statistics of u sufhciently. However, 
constructing the sample set does not require EIT measurement data; it only depends 
on the measurement setup, geometry and the prior models. The same applies to 
computing the eigenvalue decompositions, matrices B^'^^, ZJPOD g^j^^^ ^pod^ ^j^g 
approximation error statistics. Hence, all these time consuming tasks can be performed 
off-line, before starting the measurements. Once these precomputations are carried 
out, the solution of the minimization problem (|37|) is obtained with an iteration with 
reduced-order models; each iteration step necessitates only calculating the sum ((5T|) . 
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Figure 1. Schematic figure of the electrode locations and numbering, and current 
injections /f-') . 

solving the low dimensional system (P^ . (in Gauss- Newton method) computing the 
Jacobian matrix of the reduced order mapping y^'-"^, and computing the estimate 
update by solving a linear system with N unknowns. If iV << N and M << AI, 
these tasks are significantly less demanding than those needed for solving the original 
MAP-estimate ([2T1). 

4. SIMULATIONS AND EXPERIMENTS 

In this section, the EIT reconstruction with the proposed POD based reduced-order 
model is evaluated with numerical and experimental tests. The results are compared 
with the reconstructions obtained using a conventional FE approximation for the 
GEM. All the computations are carried out in Matlab environment with Dell Precision 
T7400 workstation (two quad core Intel Xeon E5420 GPUs and 32 Gb of RAM). 

4.I. Prior models and POD bases 

We modeled the target domain as a circle with the diameter 28 cm. In the model, 16 
electrodes (width 2.5 cm) were set equidistantly on the boundary. Electric currents 
with an amplitude of 1 mA were injected to the target with opposite current injection 
scheme, see Figure [T] The contact impedances ze were set to 0.01 ilcm^ for all 
electrodes. The 2D computational models described in this section were used in the 
EIT reconstructions both in numerical simulation studies and in the experimental 
studies. In the real data measurements described in Section 14.31 the target domain 
was a 3D cylinder (diameter of 28 cm, height 7 cm), the conductivity was homogeneous 
in vertical direction and the electrodes extended from bottom to top of the cylinder. 
In such a case, the 2D approximation of the EIT model is usually adequate. In all the 
inverse computations (corresponding to both conventional and reduced-order models) 
we used FE meshes consisting of 2414 elements. In the conventional FE scheme, we 
approximated the conductivity distribution in piecewise linear basis and the potential 
distribution in piecewise quadratic basis. The dimensions of the bases {<j>i{r)}iLi ^-^^d 
{ipjiT^yfLi ill the conventional FE-approximation of the GEM were = 1263 and 
M = 4939. 

In order to study the effect of the prior model to the model reduction, we 
constructed the POD based reduced-order models corresponding to two different prior 




0.5 1 1.5 2 2.5 3 



a (mS/cm) 




-0.5 0.5 1 1.5 



u{V) 

Figure 2. First row: Conductivity samples a^'\ i = 1,...,5 drawn from 
the prior TTpr,!. Second row: Potentials corresponding to the conductivity 
samples in upper figure. All potential distributions correspond to the first current 
injection, (cf. Fig. [T}. 

distributions. We denote the probability densities corresponding to the two priors by 
"■pr.ilc) and 7rpi.^2(o'). Both of the prior models were selected to be of the form of the 
proper (informative) smoothness prior [5l|2Tl[25]. In the two models, the variances 
of the nodal conductivity values ai = cr(ri) were equal, var(cri) = 0.25/iS^cm~^, i = 
1, . . . ,N. The degree of the spatial smoothness, however, differed between the two 
models; in the model Trp-[,i{cr), the cross-covariance cov(cri,crj) decreased faster with 
the distance between the coordinates r; and than in the model 7rpi..2(o')- Five 
random samples corresponding to prior models 7rpi..i(cr) and 7rpi.^2(o') are depicted on 
top rows of Figures [5] and [3J respectively. As expected, the samples corresponding to 
7rpi._2(o') are smoother than those corresponding to 7rpr^i(tT). On the bottom rows of 
Figures [2] and [3l the electric potential fields corresponding to the selected conductivity 
samples and the first current injection are illustrated. 

The POD modes corresponding to the priors TTpi- i and 7rpi..2 were constructed as 
described in Section [3] The sample sets consisted of 7500 conductivity and potential 
samples. In the selected test case, the current injections and prior distribution for 
the conductivity were rotationally symmetric. Hence, the POD basis for the potential 
was only generated corresponding to the first current injection j'"'^-'; the potential 
bases corresponding the other current injections /^*\ i = 2,..., 8 were constructed 
by rotating the first basis by 360/16 ■ {i — 1) degrees. If the current injections or 
prior distribution were not rotationally symmetric, the POD basis would have been 
generated for each current injection separately. The first five POD bases of the 
conductivity and the potential corresponding to the prior models TTpi-.i and 7rpr_2 are 
drawn in Figures |4] and [SJ respectively. As expected, the POD bases (jif^^ of the 
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Figure 3. First row: Conductivity samples a^'\ i = 1,...,5 drawn from 
the prior 7rpi._2- Second row: Potentials corresponding to the conductivity 
samples in upper figure. All potential distributions correspond to the first current 
injection. 
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Figure 4. POD bases corresponding to prior density TTpr i. First row: 
POD bases 0^*-"^, i = 1,...,5 of the conductivity. Second row: POD bases 
tpj'-'^, i = 1, . . . , 5 of the potential. 



conductivity a are spatially smoother in the case of prior TTpj. 2 than in the case of 
prior TTpr,!. The POD bases ipi of the potential u possess qualitatively similar features 
in cases of the two different prior models. 

The fractions of the variances that are captured by using A^- and il/-order 
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Figure 5. POD bases corresponding to prior density 7rpr,2- First row: 
POD bases (/i?*^^, i = 1,...,5 of the conductivity. Second row: POD bases 
^POD^ i = 1, . . . , 5 of the potential. 




N M 

Figure 6. Left; Xcti the fraction of the variance retained by the reduced-order 
representation of cr as function of TV, the dimension of the reduced order basis 
{</>?*-"-'}. Right: Xtii the fraction of the variance retained by the reduced-order 
representation of u as function of M, the dimension of the reduced order basis 
{i^JOD}. The the sohd black lines correspond to the prior model TTpr,! and the 
dashed gray lines to the prior model 7rpr,2 

approximations of a and u, respectively, are given by X(j{^) = (X^feLi -^fc)"^ Sfci ^1 

and XuiM) = {Y.k=iK)~^llf=iK- Figured] displays Xa{N) and Xu[M) for both 
prior models TTpr i and T^p-c,2- A rapid convergence of the retained variance fraction to 1 
indicates that an accurate POD approximation requires only a small number of POD 
bases. Clearly, in the case of prior 'Kp-[^2, the retained variance fraction Xa increases 
more rapidly than in the case of prior TTpi-^i. This result is intuitively appealing: when a 
random field features a very high spatial smoothness (prior model i:-pT:.2), its variations 
can be represented in a low dimensional basis. For the same reason, Xu increases more 
rapidly than Xa ^ indeed, the samples of the electric potential u are spatially much 
smoother than the samples of the conductivity a (cf. Figs. [5] and [3]). Further, the 
properties of a are reflected by the smoothness of u via the diffusion model ([1]); the 
smoother a is, the smoother is also u. In consequence, Xu increases more rapidly in 
the case of prior 7rpr,2 than in the case of prior 7rpr,i. 
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In this section, we demonstrated the effect of the prior model to the POD based 
model reduction by illustrating the samples, POD bases and retained variance fractions 
in cases of to two prior models TTp^ i((t) and 7rpr.2(cr) corresponding to different degrees 
or spatial smoothness. Below, we evaluate the POD based EIT reconstruction. The 
MAP estimates for the conductivity distributions are computed only using the prior 
model 7rpr,i(cr), that assumes higher spatial variation for the conductivity. 

4-2. Simulation results 

The simulation tests were carried out with three different target distributions: a target 
with a smooth resistive inclusion (test case 1), a target with three rectangular resistive 
inclusions with sharp boundaries (test case 2), and a target with both a resistive and 
a conductive inclusion (test case 3). The targets are illustrated in Figures [THH] (left). 

The simulated measurement data were computed using the conventional FE 
approximation of the measurement model. The FE mesh consisted of 8394 elements. 
The conductivity distribution was approximated in piecewise linear basis and the 
potential distribution in piecewise quadratic basis. The number of nodes in the 1st 
order FE mesh was 4374 and in the 2nd order mesh 17141. Gaussian distributed noise 
was added to the simulated measurements. The noise consists of two components: 
both of the components were of zero mean; the standard deviation (std) of the first 
component was 1% of the absolute value of the noiseless voltages, and the std of the 
second component was 0.1% of the difference between the maximum and minimum 
voltages. 

The EIT reconstructions with both the conventional FE approximations and the 
POD based model reduction were computed. In the former case, the reconstruction 
was computed by solving the minimization problem (|2ip . and in the latter case by 
solving (|37p. Both optimization problems were solved with the Gauss- Newton method, 
employed with a line-search [4^. The same prior model 7rpi.(fT) = 7rpi.^i((T) (see Section 
14. ip was used both in conventional and in the reduced-order reconstructions. Further, 
in both cases, the variance of the measurement noise was assumed to be known. 
The numbers of POD modes both for cr and u were selected such that 99% of the 
variances of the respective random variables were retained. This yielded reduced- 
order representations consisting of 54 basis vectors for the conductivity and 25 basis 
vectors for the potential. 

The reconstructions corresponding to both the conventional and the reduced- 
order model are depicted in Figures [7][9] (middle and right). The times required for 
each reconstruction are also indicated in the figures. In test case 1, the position of 
the inclusion was well tracked in both reconstructions. The computation time for the 
reconstruction corresponding to the conventional model was about 37 s, and for the 
reconstruction with the reduced-order model significantly less, below 200 milliseconds. 
That is, the reconstruction times were more than two orders of magnitude smaller in 
the case of the reduced-order model. In test cases 2 and 3, the reconstruction times 
were approximately the same as in test case 1. Also in these test cases, both standard 
and reduced-order reconstructions correspond well to the true targets. In all test 
cases, the POD-based reconstructions are only slighly smoother than those based on 
standard FE bases. 

In test case 1, the reliability of the POD based reconstruction was assessed 
by computing the confidence limits defined by two standard deviations (std) for 
the reconstructed conductivity. The profiles of the true and the reconstructed 
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Figure 7. Test case 1: Left: True simulated target (Ttme- Middle: Conventional 
reconstruction ctfem- Right: Reduced-order reconstruction <7poD- The 
reconstruction times are shown in the parenthesis on top of the reconstructions. 
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Figure 8. Test case 2: Left: True simulated target crtruc Middle: Conventional 
reconstruction ffpEM- Right: Reduced-order reconstruction apoo- The 
reconstruction times are shown in the parenthesis on top of the reconstructions. 
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Figure 9. Test case 3: Left: True simulated target crtruc Middle: Conventional 
reconstruction ctpem- Right: Reduced-order reconstruction o-poD- The 
reconstruction times are shown in the parenthesis on top of the reconstructions. 



Reduced-order model for EIT based on POD 



16 



OpoD 

2std limits 

-10 -5 5 10 

r(cm) 

Figure 10. Test case 1: Profiles of the true conductivity distribution (solid black 
line), the estimated conductivity irpoD (solid blue line) and the 2 std limits of 
the estimate (dashed blue lines). 

conductivities, and the 2 std limits on a cross-section of the target domain are plotted 
in Figure[Tni Here, we selected the cross-section such that it diagonally passes through 
the midpoints of the inclusion and the center of the circular target domain. The 
figure shows that the true conductivity is mostly between the 2 std limits. The figure 
also reveals that the uncertainty of the estimate is lowest near the target boundaries 
and the highest in the middle of the target. This is because the sensitivity of EIT 
measurements is highest near the target boundary where the electrodes are attached. 

To study the effects of dimensions TV and M of the POD bases {<Pi'~'^}iLi £^nd 
{iljJ'^^}jLi to reduced-order reconstructions, we computed the MAP estimates (|57| 
corresponding to test case 2 using various choices of N and M. In Figure [TTl the 
results are presented in the form of a table where increases from bottom to top, 
and M increases from left to right. The computation times are written above the 
corresponding reconstructions. The results demonstrate that especially the number 
of basis functions for the potential (M) can be set very small without ruining the 
quality of the reconstructions - even the reconstructions with M = 5 are relatively 
good, given that N is large enough {N > 45). Clearly, the number of basis functions 
required for approximation the conductivity (N) is higher than that of potential - the 
reconstructions corresponding to < 15 are severely blurred. These are expected 
results, because with the selected prior model iTp^-^a) ~ 7i'pr,i(o') the number of basis 
functions required for retaining essentially all variations of cr is A « 30, while retaining 
the variations of u requires somewhat less basis functions, cf. Figure [6l 

4-3. Experimental results 

In the experiments, a cylindrical tank was filled with tap water, and three different 
target conductivities were constructed by inserting both resistive and conductive 
objects into the tank. The first target included one plastic (resistive) bar, the second 
target two plastic bars, and the third target contained one plastic and one metallic 
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Figure 11. Reduced-order reconstructions and reconstruction times with 
different numbers of basis functions for a (TV) and u (M). The true target is 
shown in Figure |9] 



bar. All these inclusions were homogeneous in the vertical direction. Photos of all 
three targets are shown in the left column of Figure [121 The diameter of the tank was 
28 cm and the height of the water level was 7 cm. The electrode configuration in the 
tank consisted of 16 uniformly positioned, 25 mm wide boundary electrodes, and they 
extended vertically from bottom to height of the water level. In Figure [T^ the brown 
stripes on the tank wall indicate the locations of the electrodes. The measurements 
were carried out using opposite current injection scheme as in the simulations. As the 
measurement system, the SIPFIN measurement device was used [J. 

Again, the EIT reconstructions were computed both based on a conventional FE 
approximation of the CEM and with the POD reduced-order model. Same models 
were used as in the reconstructions based on simulated data in the previous section. 
The same numbers of POD basis functions were selected as in simulations studies with 
targets 1-3, i.e. iV = 54 and M = 25. 

The results are shown in Figure [T2l In the reconstructions, both corresponding to 
the conventional FE approximation of the CEM and the POD reduced-order model, 
the inclusions are clearly visible and the positions of the inclusions correspond to 
the real targets. In all test cases, the qualities of the reconstructions corresponding 
to the POD reduced-order model are comparable to those corresponding to the 



i SIPFIN is a modific ation 
|http: / /www.radic-research.de| 



of the Radic Research SIP256 instrument, 
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conventional FE based model. Again, the computation times were remarkably shorter 
when reduced-order model was used: with the conventional FE approximation, the 
computation times varied between 37903 ms and 56394 ms, while with the POD 
approximations the computations took only 169 ms in all test cases. 

Some artifacts are present in all reconstuctions. These are due to 1) modeling 
errors that were not accounted for in the observation model. Such modeling errors are 
the discretization error and the error resulting from the unknown contact impedances. 
These errors could be accounted for by the approximation error modeling [211 [32l [33] . 
2) Secondly, the selected prior model, the smoothness prior, is obviously not the best 
possible choice in the test cases. Indeed, accross the boundaries of the inclusions, the 
conductivity is discontinuous, and the conductivity contrasts are high. In such cases 
the assumption of smoothness is not a feasible one. However, as the aim of this study 
was only to demonstrate the feasibility of the POD based model reduction, we omit 
the considerations of other modeling errors and prior models from this paper. 

5. CONCLUSIONS 

In this paper, a reduced-order model for EIT has been proposed. The model is based 
on the POD representations of the electric conductivity distribution and the potential 
distribution. The model reduction is considered in the Bayesian inversion framework, 
and the POD bases are constructed based on a priori information on the conductivity. 
In the reconstructions, the errors caused by the model reduction are treated with 
the approximation error method. The reduced-order measurement model has been 
tested both simulated simulated and experimental data. The results show that by 
expressing the conductivity and the potential with reduced-order basis representation, 
it is possible to obtain feasible reconstructions with a low computational effort. In 
the selected test cases, the speedup was more than two orders of magnitude. In this 
paper, the POD based model reduction was tested with 2D examples. However, the 
proposed method can be directly implemented to computational 3D model of GEM in 
the presented form. 
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